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Abstract 



o 
u 

c I We show the existence of a series of transforms that capture several 

structures that underlie higher-dimensional partitions. These transforms 
lead to a sequence of triangles whose entries are given combinatorial in- 
terpretations as the number of particular types of skew Ferrers diagrams. 
The end result of our analysis is the existence of a triangle, that we denote 
by F, which implies that the data needed to compute the number of parti- 
tions of a given positive integer is reduced by a factor of half. The number 
of spanning rooted forests appears intriguingly in a family of entries in the 
F. Using modifications of an algorithm due to Bratley-McKay, we are able 
to directly enumerate entries in some of the triangles. As a result, we have 
been able to compute numbers of partitions of positive integers < 25 in 
any dimension. 
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1 Introduction 

An unrestricted <i-dimensional partition of n is a collection of n points (nodes) 
in Z + +1 satisfying the following property: if the collection contains a node a = 
(ai,a2, ■ ■ ■ ,a,d+i), then all nodes x = (xi,x 2 , . . . ,Xa+i) with < Xi < a» V i = 
1, . . . , d + 1 also belong to it [UI2]. Let p<j(n) denote the number of distinct such 
partitions. Denote by Pd{q), the generating function of unrestricted <i-dimensional 
partitions: (pd(0) = 1) 



oo 



P d (q) = }^ p d (n) q n . (1.1) 

n=0 

There exist explicit formulae for the generating functions for d — 1 and d = 2 
due to Euler and MacMahon respectively [3]. However, no such formulae exist for 
d > 2 as an inspired guess of MacMahon was subsequently proven to be false p] . 
It appears that there is no simple formula and one has to take recourse to brute 
force enumeration. Given that asymptotically one has [IH6] 

logp d (n) ~ n d ' d+1 , (1.2) 

it is easy to see that the numbers grow exponentially fast and naive enumeration 
is not the way to go. 

The first serious attempt at direct enumeration of partitions in any dimension 
is due to Atkin et. al. [lj based on an algorithm due to Bratley and McKay 
[7]. Knuth provided another algorithm that enumerates numbers of topological 
sequences which can be used, in principle, to generate numbers of partitions in 
any dimension [8]. Both algorithms are highly recursive and easily implemented 
on a computer. 

This paper attempts to find structures in the enumeration of partitions and 
come up with refinements in their enumeration. Such refinements when cleverly 
combined with computer-based enumeration should in principle enable one to 
enumerate partitions of integers below some maximum value in any dimension. 
The maximum value turns out be 25 in our case though we believe that, with 
some effort, this number can be pushed to around 30. 

Our refinements begin with the result of of Atkin et. al. who showed that the 
binomial transform of Pd(n) leads to a lower-triangular matrix that we denote by 
A = {a n ,r)- 

d+l , , 1 \ 
Pd{n) = X) ( J dn,r ■ (E3D 

This transform implies that in order to compute partitions of a positive integer 
n in any dimension, we need to only compute (n — 1) numbers that make up 
a particular row of the triangle A. We show the existence another triangular 
matrix, that we denote by F = (f n ,x)i as a transform of the matrix A with fewer 
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entries. 

r m 

Q>m+r+l,r = / / 

x=0 p=x 

Our result is that we need only [(n— l)/2] independent numbers i.e., roughly half 
of the initial estimate to determine partitions of n in any dimension. We illustrate 
the gain by explicitly displaying the first eleven rows of the A and F-matrices. 
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The F-triangle is, in a sense, the end-point of a sequence of transforms and 
triangles that we introduce. We also provide combinatorial interpretations for the 
various triangles that appear as a result of these transforms. This enables use to 
modify the Bratley-McKay(BM) algorithm to directly enumerate the matrix A 
that we mentioned earlier and a second triangle, C that we define in the sequel. As 
we discuss in the appendix, similar refinements can be carried out for partitions 
restricted in a box. 



1.1 Summary of results 

1. Given a partition in any dimension, we have introduced two new attributes: 
its intrinsic dimension (i.d.) - see definition 12.11 and its reduced dimension 
(r.d.) - see definition 12.71 

2. These two attributes lead to two new triangles, the A and C-matrices (see 
Eq. f!2.4p and ( 12. lip ) whose entries admit combinatorial interpretations. 
We propose a further refinement in the form of two other triangles, the D 
matrix(see Eq. (" 12. 15[) ) . 

3. We show that the C/D triangles are the first in a series of transforms, the 
end-point of which leads to a triangle F (see Eq. ( I2.2(jp ). The n-th row 
of this matrix has only [(n — l)/2] entries (where [x] is the integral part of 
x) and these entries determine the partitions of n in any dimension. This 
constitutes the main result of this paper. 

4. We see an intriguing relationship between the numbers of spanning rooted 
forests on m vertices and a components and a family of entries in the F- 
matrix. This is Proposition 12.161 

5. We conjecture the existence of two other triangles, the a- and the /3-matrices 
with integer entries. 



6. We prove a conjecture of Hanna on the existence of a triangle that deter- 
mines all higher-dimensional partitions. 

7. We propose a modification to an algorithm of Bratley and McKay that 
enables us to directly compute the A and C matrices. We compute the first 
25 rows of the F-matrix thereby obtaining partitions in all dimensions for 
integers < 25. 

8. Tables [THE] provide the numerical results that we have obtained. 

2 Structures in higher-dimensional partitions 

2.1 Ferrers diagrams and permutation symmetry 

A Ferrers diagram represents the partition as a (d + l)-dimensional arrangement 
of nodes. For instance, the following one-dimensional partition of 4 

or ( o o ? 2 ) in compressed form 



oy ' \oj ' \ij ' ^2 

is represented by the following two-dimensional Ferrers diagram or as a Young 
diagram where we replace the nodes by squares(more generally, hypercubes). 

%2 
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There is a natural action of Sd+i on the (d + l)-dimensional Ferrers diagram - 
this corresponds to permuting the (d+1) coordinates. For one-dimensional parti- 
tions, this is referred to as conjugation. The symmetry group of a d- dimensional 
partition is the largest sub-group of Sd+i that acts trivially on the corresponding 
Ferrers diagram. 

2.2 The intrinsic dimension 

Typically, one is interested in the asymptotic behavior of Pd(n) for large number 
of nodes n while keeping the dimension d fixed. However, one may ask about what 
happens to Pd{n) if we keep the number of nodes, i.e., n, fixed and keep increasing 
d. It is easy to see that when d > n + 1, all the nodes of the Ferrers diagram 
necessarily lie in some r-dimensional hyperplane with r < d. This motivates the 
following definition (implicitly present in Atkin et. al. [1]). 

Definition 2.1 Given a Ferrers diagram, let it be contained in a r-dimensional 
hyperplane but not in any (r — 1) - dimensional hyperplane. The intrinsic dimen- 
sion(i.d.) of the Ferrers diagram is defined to be r. 



Note that such a r-dimensional hyperplane is given by setting (d + 1 — r) 
coordinates to zero. Any permutation of the (d + 1 — r) coordinates (that are 
set to zero to obtain the hyperplane containing the nodes) does not change the 
Ferrers diagram. It is thus easy to see that the symmetry of a Ferrers diagram in 
(d + l)-dimensions of i.d. r is necessarily of the form H x Sd+i~ r where H C S r . 
We shall (somewhat loosely) call H, the symmetry of the Ferrers diagram. 

Let two d- dimensional partitions be equivalent if their Ferrers diagram are 
related by an S d +i action. It is easy to see that all rf-dimensional partitions 
belonging to such an equivalence class have the same intrinsic dimension. Further, 
given a d + 1-dimensional Ferrers diagram with symmetry H and i.d. r, the 
number of Ferrers diagrams in its equivalence class is given by the order of the 
coset S d+1 /(H x S d+1 _ r ) i.e., 

(d+l)\ /d+T 

- ' 'X 



(d + 1 — r)\ x ord(iJ) \ r / ord(H) 



Definition 2.2 A Ferrers diagram is said to be strict when its intrinsic dimen- 
sion equals its dimension. 

Given a d + 1-dimensional Ferrers diagram of i.d. r, it is useful to drop the 
(d + 1 — r) dimensions that are orthogonal to the hyperplane containing the 
nodes thus obtaining a strict FD. The symmetry of the strict Ferrers diagram is 
now H C S r . 

Definition 2.3 A generalized Ferrers diagram (gFD) refers to the equivalence 
class of strict Ferrers diagrams obtained by the action of S r on a given strict 
Ferrers diagram of i.d. r. 

The number of strict FD's in a gFD of i.d. r and symmetry group H is or .^ H ) ■ 

Definition 2.4 The weight of a gFD of i.d. r and symmetry H C S r is defined 
to be r, ! m . 

ora(H) 

Since H C S r , Lagrange's theorem implies that the weight, L H \ , is a positive 
non-zero integer. Note that the weight is independent of the dimension of the 
Ferrers diagram and is the same for all elements in an equivalence class. Thus, 
to an equivalence class of a given Ferrers diagram, we associate three numbers: 
the number of nodes n, the i.d. r, and the weight, w. An important observation 
is that there exist no Ferrers diagram with n nodes and i.d. r > n - this follows 
from noting that one needs at least r + 1 nodes to create a Ferrers diagram of 



i.d. r. We see that the number of d- dimensional partitions is thus given by 



Pd{n) 



r=0 V J X\-(n,r) 

r=0 V J [A]h(n,r) 

n-1 
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r=0 



d+1 



r 



a 
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(2.1) 
(2.2) 
(2.3) 



where the second line defines a n , r as the sum over all strict FD's with n nodes and 
i.d. r. In the second line, the sum over [A] indicates that we sum over equivalence 
classes of strict Ferrers diagrams (gFD). Note that a n>r has no dependence on d 
and counts the numbers of strict Ferrers diagrams with n nodes and i.d. r. We 
shall provide a second, and more useful, combinatorial description of a nr later. 



2.3 The first transform 

We extend a nr into a lower-triangular matrix, that we denote by A, by setting 
a nr = when r > n. Thus, we obtain the matrix A = (a nr ) for n — 1,2,... and 
r = 0, 1, 2, . . .. With this definition, we an rewrite the above equation as 



(2.4) 



To our knowledge, the above observation first appeared in a paper by Atkin et. al. 
PP. Thus the Pd{n), for a fixed value of n, corresponds to the Binomial Transform 
of the n-th row of the matrix A. It is easy to see that a n> o = £ n ,i- The lower 
triangular nature of A implies that only (n — 1) numbers, (a n ,i, On,2, • • • , Q>n,(n-i)) 
determine pd{n) for any d. The matrix A appears in the OEIS as sequence number 
A119271 [9]. The inverse Binomial transform is given by 





(2.5) 



with po{n) = 1. Of course, a nr = when n < r+ 1 reflecting the lower-triangular 
nature of the matrix. Suppose we know all partitions of n max up to d max . This 
determines the first n max rows and (d max + 1) columns of the matrix A. 

For low values of n, we can explicitly compute the entries in the A-matrix by 



listing the gFD's and working out their weights as we do below. 



P*(2) = (T) 

Pd(3) = ft 1 ) 
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(2.6) 

(2.7) 
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= (r)+3(r)+(T)- 

The first few rows of the A-matrix are as follows (see also [TD] 
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Definition 2.5 Consider a pair of FD's (A,/i) such that \x C A. Then, a skew 
Ferrers diagram is the set of nodes X\fi. 

One can think of the entries in the A-matrix as counting skew Ferrers diagrams 
obtained by deleting the node at the origin (0,0,..., 0) T that is contained in any 
Ferrers diagram. Then, a n , r is the number of strict FD's of dimension r obtained 
by adding (n — 1) nodes to the node at the origin. One sets a\$ = 1. 



2.4 A combinatorial interpretation 

We will now provide another combinatorial interpretation for the numbers a nr 
that make up the lower-triangular matrix A. We begin with the observation that 
a r+ i tr = 1 - this follows because there is a unique FD of i.d. r containing r + 1 
nodes. The coordinates are given in the following r x (r + 1) matrixLj 



10 
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(2.9) 
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1 Recall that each column is the coordinate of a node and thus there are (r + 1) columns and 
r rows. 



This FD has maximal symmetry S r and weight 1. 

Remark: Every FD with intrinsic dimension r necessarily contains fi r . This 
implies that an FD with n nodes and i.d. r can be obtained by adding m = 
n — r — 1 additional nodes to fi r . This leads to the following combinatorial 
interpretation for a m+r+ i jr .. 

Proposition 2.6 a m+r+ i yr is the number of strict Ferrers diagrams with i.d. r 
obtained by adding m nodes to the standard Ferrers diagram, \x r . 

Let A be an FD that contributes to a n , r . Its symmetry group H C S r - this 
implies that there will be r\/odd(H) = wt(X) distinct FD's obtained from it by 
the action of S r . It is easy to see that the process of adding m nodes to fi r will 
generate the same number of FD's that belong to the equivalence class (gFD) [A]. 
So far we have completely determined the first 25 rows of the A-matrix (see 
Table [2]). The entries have been determined by combining several methods: (i) 
taking the inverse Binomial transform of known numbers for higher-dimensional 
partitions, (ii) by direct enumeration using the combinatorial interpretation and 
(iii) by determining another triangle, the C-matrix, that we introduce later. It is 
important to note that the numbers, when available, from the different methods 
agree. Further, none of the conjectural formulae are used in determining the 
entries. 

2.5 The second transform 

Definition 2.7 Let A be an FD of i.d. r and consider the skew FD A \ /i r . Let 
the nodes of the skew FD be contained in a x-dimensional hyperplane (obtained by 
setting r — x coordinates to zero) but not in any (x — 1) -dimensional hyperplane. 
The reduced dimension (r.d.) of the FD A is said to be x. 

Clearly the reduced dimension of an FD is always less than or equal to its intrinsic 
dimension. The symmetry of a FD with i.d. r and r.d. x is necessarily of the 
form H x S r - X C S r . Then, one has 



x=0 



m+r+l,r / I I Cm, a: > v ^v 



where the binomial term ( r ) takes into account the situation with maximal sym- 
metry and co,o = 1 and c m $ = co, m = for m > 0. 

1. The coefficients c m)X are clearly independent of the the i.d. (r) as they are 
related to the skew FD's with m nodes and r.d. x. 

2. We say that a skew FD is strict if its dimension and r.d. are the same. 



3. Let us denote the equivalence class of strict skew Ferrers diagrams, A \ fj, x , 
under the S x action as an sFD. All skew FD's in an sFD will have identical 
reduced and intrinsic dimensions. Thus, given such a skew Ferrers diagram 
with symmetry H C S x , its equivalence class will contain ov ^ H \ distinct 
skew Ferrrers diagrams. 

4. The c m)X are non-negative integers since they count the number of strict 
skew FD's with m nodes and r.d. x. 

5. For fixed m, one can see that the maximum value of r.d. with m nodes is 
2m. This enables us to convert the above equation into a second binomial 
transform 



(2.11) 



where we extend c m>x into a triangle, C = (c mtX ), by setting c m>x = for 
x > 2m. We usually do not write out the zeroth row and column of the 
C-matrix. 

6. For fixed m, we can consider a m+r+ i^ r as a function of r. The function 
9m{f) '■= 2m!! a m+r+ i )r is a polynomial of degree 2m, conjecturally with 
integer coefficients, in the variable r and g m (0) = for m > 0. 

7. We have directly determined eleven rows (m G [0, 10]) of the C-matrix (see 
Table |3]) . The first few rows of the C-matrix are: 
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It is easy to see that there is only one sFD with m nodes and r.d. 2m. In the 
picture below, the m nodes of the sFD are indicated by open circles. The filled 
circles indicate the nodes of /x 2m that must be added to the sFD to obtain an FD. 




X \° 

^2m-l 



(2.12) 



The symmetry of the skew FD is (S m x Z™) and thus c m)2m is the dimension of 
the coset i.e., 

_ dim(ff 2m ) _ 2m! _ 

m ' 2m ord{S m x Z?)) 2m!! l j " ' 

Definition 2.8 A skew FD of i.d. r is said to be reducible if a proper subset of 
its nodes are contained in a d-dimensional hyperplane (obtained by setting r — d 
coordinates to zero) with d < r and the nodes not in the proper subset lie in 
the orthogonal (r — d)- dimensional hyperplane (obtained by setting the other d 
coordinates to zero). 

Definition 2.9 We say that an FD, X, of i.d. r is reducible if the skew FD, 
X \ (x r is reducible. 

Thus a reducible sFD has multiple components consisting of non-intersecting 
proper subsets of its nodes lying in mutually orthogonal hyperplanes. Thus the 
sFD given in Eq. f 1 2 . 1 2 j) is reducible with m components each of which is isomor- 
phic to the irreducible sFD a 2 defined as follows: 



a 2 = A 



(2.13) 



We can thus write the sFD (j2.12p as o 2 x a 2 X • • • x cr 2 = a™. 

Similarly, one has two distinct sFD's with x = 2m — 1 and the two sFD's are 
reducible containing a 2 (for some suitable value of n) as one of the components 
and the other component are the following two irreducible sFD's that contribute 
to c\ i and c 2 3 respectively. 



_ a 3 = xT° = ( 11 

01= • » O > Xn^ ^ 01 



(a) (b) 



(2.14) 



where we have called the second sFD 03 - it has two nodes and has r.d. 3. In other 
words, c 2m:2m -i has contributions from two sFD's - one of the form a 2 x o\ 
and the other of the form 2 ~ x a 3- Studying the symmetries of these two 
sFD's with r.d. (2m — 1), one obtains 

(2m -1)! (2m - 1)! (2m -1)! 

Cm,2m-1 = -77, ^TTT + 7777 7TTT = m X 



(2m -2)!! 2(2m-4)!! (2m - 2)!! ' 

Clearly, such a diagrammatic method will enable one to write further formulae 
(we will provide a few more in an appendix) for c mjX . However, it can get tricky to 
find all possible diagrams. Keeping this in mind, we make the following definition. 



Definition 2.10 The density, p, of a sFD with m nodes and r.d. x is p = m/x. 

The density of a sFD is always greater than or equal to | since c m>x = when 
x > 2m. 

Proposition 2.11 When its density is in the range (|, |), an sFD with m nodes 
and r.d. x is necessarily reducible and one of its components is the sFD, (cr 2 ) n , 
for some n > n min = 2x — 3m. 

The proof follows from Proposition 12.151 that we prove later. When p < 2/3, the 

proposition implies it is impossible to construct an sFD that does not contain o-i 

as a component. The first new sFD, 03, appears at p — |. The minimum value 

of n is fixed by the condition that the density of the sFD goes past or equals | 

after deleting the nodes that appear in (02)™ i.e., it is smallest value of n such 

that 

m — n 2 

> - ==>- n >2x — 3m . 

x — 2n 3 

2.6 The third transform 

Proposition 12. Ill suggests that in counting the skew FD's that contribute to c m , x , 
we can remove components isomorphic to 02 in reducible skew FD's and only 
count skew FD's that do not contain any 2 components. This motivates the 
next transform where we introduce a new triangle D = (d mtX ). 



2^ (Oo,\\\( V -0n,\\ m -i 
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(2.15) 



with d fi = 1, d m fi = c?o,m = for m > and y min = 2x — 3m. The pre- factor 
in the transform is determined by the order of the symmetry of a\ which is 

2yy\ = (2y)\\. 

1. d m ^ x counts the number of skew FD's with m nodes and r.d. x not containing 
o~2 as its components. Thus it is positive definite. 

2. Proposition 12.111 implies that d m ^ x = when m/x > 2/3. This is stronger 
than the condition m/x > 1/2 implied by the property of the C-matrix. 

3. It is useful to rewrite the transform as follows: 

E(2m — z)\ 
(2m - 2u) U(2y - z)\ 

In this form, one sees that completely determining row z of the D-matrix 
leads to a nice compact formula for c m ^ m -z- The D-matrix clearly contains 
fewer terms than the C-matrix since d mjX = when p < 2/3. 

10 



4. To illustrate the transform, consider c m ^m-i which we have already com- 
puted. One sees that 
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c 2m , 2m _i - 2J (2m ^ 2j/)!!(2 y j/ - 1)! rfj/,2y - 1 
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(2.17) 



It is easy to see that d\ t \ = 1 as there is precisely one sFD ((a) in Eq. (12.141) ) 
and g?2,3 = 3 as there are three inequivalent diagrams under the action of 
5*3 on the sFD, cr 3 . 

5. When p = 2/3, there is only one sFD, a™, that contributes to rf 2 m,3m- This 
implies that 

3 "' ! (2.18) 
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2.7 The final transform 

The main advantage of the .D-matrix is that it contains fewer terms than the 
C-matrix. Using it, we have arrived at formulae for c m ^ m -z for z = 2,3,4,5 
analogous to the one in Eq. ( I2.17P that can be obtained, in principle, from the 
table that gives the .D-matrix. Can we do better? We saw that as the density 
increased from 1/2 to 2/3, only one irreducible diagram appears. At p — |, two 
new sFD's appear. They are 



C"4a 
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1 
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10 
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(2.19) 



In fact, one can define another transform that removes reducible components of 
type 0"3 from sFD;s that contribute to the D-matrix for p G (2/3, 3/4). The next 
proposition will enable to do this and a lot more by removing a whole family of 
reducible components that necessarily appear when in sFD's with p < 1. 

Definition 2.12 Let T> = U r V r , where T> r denotes the set of strict Ferrers dia- 
grams of dimension r consisting only of nodes of the form (1, 1,0,..., 0) T or its 
S r images in addition to the nodes present in p r . 
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We say, somewhat loosely, that a strict skew FD, a of r.d. x is in T> if the FD 

fi x U a G T>. One can show that a 2 , <j 3 and cr^/t, are the only irreducible strict 

skew Ferrers diagrams at dimensions 2, 3 and 4 respectively that appear in T>. 

Let e myr denote the number of Ferrers diagrams in T> obtained by adding m 

f( x )\ 
nodes to \x r . It is easy to see that e m>x — I J as there are (*) possible nodes 

from which we need to choose m nodes. We define a new transform that removes 
reducible components that are in T>. 



r m 
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xj \m — p 



Jp+x+i, x , (2.20) 



x=l p= 

where in the second line we use the explicit formula for e m ^ x and / 10 = 1, f n ,o — 
/i,n-i — for n > 1. In the first line, a typical term in the summation on the 
right hand side consists of reducible strict FD's with the component in T> having 
i.d. r — x and (m — p) nodes added to fi r -x and the other component consisting 
of an strict FD with no reducible component in V, i.d. and r.d. x and p nodes 
added to jj, x - their number is counted by f p + x +i, x - The binomial factor ( r ) is the 
number of ways one can choose x dimensions occupied by the FD's contributing 
f° fp+x+i,x- The above formula defines a new triangle F = {f n ,r)- The entry 
fr+m+i,r is the the number of strict FD's of i.d. r obtained by adding m nodes 
to fi r and does not contain any reducible components that are in T>. Such an FD 
must necessarily have r.d. also equal to r, else it will necessarily have a reducible 
component isomorphic to fi r -x if its r.d. is x. 

It is easy to see that f r +i,r = 0. The only contribution to a r +i,r is the unique 
FD /x r which is T>. Similarly, f r +2,r = when r > 1 as the only contribution to 
a r+2,r is of the form a\ x a r 2 ^ 1 . One also has fs t i = 1 with o~\ being the unique 
FD contributing to it. The next proposition shows the advantage of defining the 
F-matrix. 

Proposition 2.13 f m + r +i,r — when r > m. 

Proof: Let A be an FD of i.d. r with m + r + 1 nodes that contributes to 
fm+r+i,r- Consider the skew FD, A \ /J, r - it has m nodes. It must be a strict 
skew FD else it has a irreducible component isomorphic to fi x for some x < r. 
Thus, the proposition implies that there are no strict skew FD's with density 
p = m/r < 1. 

We can also assume that the skew FD is irreducible - if it is reducible, it 
must necessarily have at least one irreducible component with density < 1 and 
we can focus on (proving the non-existence) such irreducible components. Our 
goal is thus reduced to proving that there are no irreducible strict skew FD's with 
density < 1. 
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Definition 2.14 Let us call the nodes obtained by all permutations of the coor- 
dinates of the node (1, 1,0,..., 0) T as nodes of type 1. Similarly, call the nodes 
obtained by permuting coordinates of (2, 0, ... , 0) T as type 2. Nodes of type 3 are 
nodes that are not of type 1 or 2. 

Examples of type 3 nodes include (1,1,1,0,..., 0) T and (3, 0, . . . , 0). Such nodes 
cannot be added to the FD /i r without including supporting nodes of type 1 and 
2. The addition of nodes of type 3, when possible, never increase the r.d. of 
an FD thus increasing the density. Thus, given a FD A (of i.d. r and r.d. r) 
containing type 3 nodes, we can form a new FD A' with the same r.d. but lower 
density. Further, if A \ [i r is irreducible, A' \ fi r is also irreducible. The skew FD 
A' \ fi r thus consists of nodes of type 1 and type 2. If it consists of only nodes of 
type 1, then A' G T>. Thus, we need to only consider irreducible strict skew FD's 
containing at least one node of type 2. 

For the rest of the discussion, let A' be an FD such that A'\/i r is an irreducible 
strict skew FD containing only nodes of type 1 and at least one node of type 2. 
It is easy to see that removing of node of type 2 does not affect the irreducibility 
of the skew FD. Further, it does not reduce the r.d. as the only way a type 2 
node can reduce the r.d. of a skew FD is when it appears as a part of a reducible 
component isomorphic to <j\. Thus, we can delete all type 2 nodes to obtain a 
new FD A" that is irreducible and contains only type 1 nodes. Again, it is easy 
to see that p(A") < p(A'). Further A" G V. Thus, one has the sequence 

p{\") < p(A') < p(A) . (2.21) 

Let A' \ fi r have (r — l)-nodes so that its density is just below one and contain 
z nodes of type 2. Then, A" \ fi r will have (r — 1 — z) nodes and be irreducible. 
The next proposition shows that such a A" does not exist. Hence, there exists no 
FD A" and hence no FD A' with density < 1. □ 

Proposition 2.15 The only strict FD's in T> of i.d. r such that the skew FD 
A \ fi r is strict and irreducible with density less than 1 necessarily have p = r — k 

Proof: Let us assume that X\ fi r has (r — 2) nodes and is irreducible. Let us 
try to construct such a strict skew FD and we will see that there are not enough 
nodes. Start by putting the first type 1 node in the x 1 x 2 plane. The irreducibility 
condition implies that the second node must be either in the x 1 x a or x 2 x a plane 
where a is not 1 or 2. The key point is that the additional node must contain 
one of the used up coordinates, x l or x 2 in this case and a new coordinate so that 
irreducibility is maintained. Clearly, such a process needs (r — 1) nodes to get an 
irreducible skew FD A \ fi r with r.d. r. This is impossible. Hence, there exists 
no irreducible skew FD A with density r ^. it is easy to extend the argument to 
exclude even lower densities. Thus, the only possibility that is not ruled out is 
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to have strict skew FD's with (r — 1) nodes with r.d. r - these have density — -. 

Remark: o~2, o"3 and a± a /b are the only irreducible strict skew FD's with r.d. 
2, 3,4 respectively. 

2.7.1 Properties of the F-matrix 

1. The most important property is the one implied by Proposition 12. 131 which 
says that the F matrix is lower triangular with f nr = when r < [(n— 1)/2]. 
For fixed value of n, the F-matrix has far fewer terms (roughly half) than 
the corresponding row in the A-matrix. We have determined the first 25 
rows of the F-matrix (see Table HJ). 

2. It turns out that there are other transforms that also lead to matrices with 
fewer entries like the F-matrix. See for instance, the box transform that we 
consider in the appendix. However, their relationship to A is not as simple 
as Eq. (j2.20p . The simplicity of Eq. (J2.20P is what picks out the F- matrix 
as special. 

3. We can also use this idea to refine the counting problem associated with 
the C-matrix. Let C v = (c^ x ) denote the contribution to the C-matrix 
that arise from FD's that are in T>. Since the set T> r is invariant under S r , 
it is easy to see that C v is given by the transform 

1m 



m I * — ' \x 



EC- ( 2 - 22 ) 



Then, we can define C = (c m)1 ) by removing contributions that arise from 
reducible parts that are isomorphic to contributions to C v . Then, one has 



x— 1 m— 1 



x — -L in—L / \ 
c m,x = c m,x + c m ,x ' / , / 4 I J C m-p,x-y C P,V • [Z.Zo) 

y =\ p =i \y/ 

Given a strict skew FD that contributes to c mtX , it is easy to see that there 
is a unique FD obtained by adding nodes in /i x to the skew FD. Further, 
this FD must contribute to the entry f m+x+1 >x in F. Since the converse also 
holds i.e, given a strict FD of i.d. x that contributes to the F-matrix, the 
skew FD obtained by deleting nodes in \i x gives a skew FD that contributes 
to C. Thus, one has 

Cm,x Jm+x+l,x ■ \^"^V 



We observe numerically that / 2m+ i ;m = (m + l) m ~ 2 for m = 0, 1, 2, ... , 12. 
We will show that it holds for all m. These numbers appear in the sequence 
numbered A000272 in the OEIS [9]. The next proposition presents a further 
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refinement. We need a few definitions which we briefly state. A graph, consisting 
of vertices and undirected edges, with no cycles is called an acyclic graph or a 
forest. A forest may consist of disconnected components and is called a tree if it 
has only one connected component. A rooted tree is one with a marked/special 
vertex (called the root) while a rooted forest is one in which every component 
is rooted. A spanning forest is any subgraph that is both a forest (contains no 
cycles) and spanning (includes every vertex) [TT|[T2]. 



Proposition 2.16 Let a be the number of nodes of type 2 contained in an FD 
that contributes to f2m+i,m- Let /2m+i,m(ct) denote the total number of such Fer- 
rers diagrams. Then, /2 m +i, m («) is the number of spanning rooted forests on m 
vertices and a components. It follows from a result due to Cayley on the numbers 
of spanning rooted forests that |I] 



/ 2m+ i, m («) = (™_ J m m ~ a . (2.25) 

Proof: We will provide a bijective map relating FD's that contribute to /2 m +i, m («) 
to spanning rooted forests on m vertices and a components. There is a natural 
action of S m on both sides - on the FD side, it corresponds to permuting the 
m coordinates and on the rooted forest side, it corresponds to relabeling the m 
nodes. We identify these two groups. 

Given a skew FD that contributes to /2 m +i,m(«), we can construct a graph 
with m vertices labeled from (1, . . . , m) as follows. The type 2 nodes become 
root vertices carrying the label of the non-vanishing coordinate. Thus if a type 
2 node has non- vanishing j-th coordinate, assign it the label j. Add (m — a) 
vertices and label them with the unused labels. Every type 1 vertex has two 
non-vanishing coordinates, say the j-th and fe-th coordinates. Assign an edge 
that connects vertex j to vertex k. Repeat for all type 1 nodes. In this process, 
there are as many components as there are type 2 nodes. Thus the graph is a 
spanning rooted forest on m vertices and a components. The following example 
illustrates the map for m = 4 and a = 1. The root vertex is shown by a filled 
circle. 



2 110 

10 

11 

1 



i > O 2 Q 3 



D4 



To prove the converse statement, given a spanning rooted forest with m 
vertices and a components, we need to construct an FD that contributes to 
/2m+i,m(«)- This is easy to do . Pick the root vertices and assign them to type 
2-nodes whose non- vanishing coordinate decided by the label of the vertex. Next 
assign to all edges a type 1 node that has non-vanishing coordinates at precisely 
the locations decided by the labels of the vertices it connects. We thus recover 
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the skew FD. □ 

An example: We know that /s ; 2 = 3. The three skew FD's are: 



0-7 



2 en . (2 i\ . ( 1 0" 

2 J 1 loiJ 1 U2, 



where the nodes are listed by the ordering: (ai, 02) > (&i, ^2) if 01 > fri or 01 = 61 
and a 2 > b 2 . Note that there are two equivalence classes of skew FD's. Under S 2 
action as the second and third skew FD's get mapped to each other. 



12 12 12 
• • : • o ; o • 



Remark: Given a skew FD, it is possible to uniquely label the nodes by ordering 
them by a choice of ordering as illustrated above. 

3 Other triangles 

3.1 New triangles 

So far, we have considered transforms that lead to new triangles (A/C/D/F) all 
of which have positive definite entries since we they all count numbers of skew 
Ferrers diagrams. We will now provide two other transforms that are partly 
conjectural and lead to triangles that are not positive definite - we denote the 
entries with Greek letters to remind us of this. We begin by expanding the entries 
in the A-matrix as follows. Let 



2m 



r 

2 = 



2m— 2 



"ra+r+l,r / (%m,z „ n j V / 



with a m fi = 1 for m > and a m ^ m = for m > 0. The above transform provides 
the entries for another triangular matrix, a myZ , that we call the a-triangle by 
setting a m , z = for z > 2m. One can explicitly relate the a m ^ z to the entries in 
the C-matrix using Stirling numbers of the first kind. Thus the above formula is 
not conjectural. However the following is conjectural: 

Conjecture 3.1 The entries of the a-triangle, i.e., ot m)Z , are all integers. 

This is true for the first ten rows and appears to hold for the first eleven rows 
which have been determined using conjectures. 

The second conjecture introduces a new triangle, that we call the /3-triangle, 
and its associated transform. It has been determined experimentally and verified 
to hold to the extent possible. 



16 



Conjecture 3.2 The a-matrix admits the following decomposition. 

a m ,z= J2 \ z ™ v ) &** ' (3 ' 2) 

y=0 V *' 

wit/i /3o,o — 1 an< ^ /?2y,j/ = /or all y > 0. 

By setting /3 2)2/ = for y > \_z/2\ , this becomes the binomial transform 

s/=o v y/ 
The inverse transform is 

/^ = E(-i) m WV) a ^- (3 - 4) 

We now state a conjecture of Meeussen that fixes one of the coefficients. 
Conjecture 3.3 (Meeussen) 

where H n (x) is the n-th Hermite polynomial. 

Recall that the a-matrix has 2m non-zero entries in the m-th row. The /3-matrix 
has fewer terms, roughly half the entries in the a-matrix. We were able to 
determine eleven rows of the a and C-matrices using the /3-matrix of which 10 
were verified through other means. This was our main motivation in searching 
for and find the combinatorial problem that eventually lead to the F-matrix. 

3.2 The B-triangle 

We now construct another lower triangular matrix B = (b n , r ) with n = 1,2,... 
and r = 0, 1, 2, . . . and b n $ = 1. 

n—l / r\ n—l 

Pd{n) = T( ) b n , r = 1 + Vri b n , r . (3.5) 



\ / r=1 \ / 



r=0 

The matrix B appear in the OEIS as sequence number A096806. Using Pascal's 
identity 

' d+A fA + ( M, (3.6) 



nr* f \ nr> f \ ty \ 

we can relate the matrix B to A. Thus, one has the relation 

One can easily show that fe ni „_i = 1 using the above formula and known properties 
of the matrix A. The first six rows of B have been determined explicitly, for 
instance, in Andrews' book on Partitions [2]. It is easy to check that the above 
relation holds for all six rows. 
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3.3 Hanna's triangle 

Conjecture 3.4 (Hanna) There exists a lower-triangular matrix T = (r^) (with 
i,j = 0, 1, 2, . . .) with integral entries and ones on its diagonal such that 

n 

Pd (n) = ^(T d ) nj . 

3=0 

In other words, the sum of the the n-th row of the d-th power of T give the 
(^-dimensional partition of n. This matrix appears in the OEIS as sequence 
A096651. Since pd(0) = 1, we can set To,o = 1 and Tjfl = for j > 0. For 
the rest of the discussion, we will consider n > and can delete the zeroth row 
and column of the T-matrix as they no longer play a role. We shall however use 
the same symbol T to denote the modified matrix as it is easy to reconstruct the 
original T matrix by adding back the zeroth row and column. We shall prove the 
existence as well as the integrality of the matrix T by constructing an explicit 
map that relates T to the matrix B (and hence A) that we considered in the 
previous section. 
Proof: For n > 1, the Hanna conjecture can written as 

n 

Pd{n) = / \J- )n,j — 2_^ T n,x 1 T x 1 ,x 2 ' ' T x d _ 1 x d , (3.8J 

3=1 xi— x d 

where n > x\ > x% > • • • > Xa > 1. It obviously holds for n = 1 since 
rii = 1. Using the fact that T has ones in its diagonal, we can simplify the above 
expression to 

P d (n) = 1 + Y^ ( J Yl r «^i r *l.«a ■ ■ • 7ir-l,Xr • ( 3 ' 9 ) 

r=l ^ ' x\---x r 

with sum now running over all sequences of r positive non-zero integers (a?i, ■ ■ • , x r ) 
such that Xo = n > X\ > x-i > ■ ■ ■ > x r > 1. The combinatorial factor expresses 
the number of ways in which diagonal elements are chosen. Comparing the above 
equation with Eq. (13. 5p implies the (potential) identity for n > 1 and r > 1. 

/ , Tn,xiTxi,X2 ' ' '1~x r -i_,x r "n,r i [Q.i-O) 

x\---x r 

with n > X\ > x-i > ■ ■ ■ > x r > 1. Let us assume that this relation holds for 
n < m) (for some m > 1) and that we have determined (m — 1) rows of T. Then, 
we can rewrite the above equation as 

J^ T m ,A,(r-i) = 6m,r for m > r > 1 . (3.11) 

l<x<m 



The above (m — 1) equations are linear equations in (m — 1) unknowns: 
(r m> i, . . . ,r mm _i) - these are the undetermined entries in the m-th row of T. 
Hence, they have a solution if the matrix (constructed using 6 Xl ( r _i)) is invertible. 
The matrix is upper triangular with ones in its diagonal. Hence it is has deter- 
minant one and hence is invertible. This enables us to recursively determine all 
the entries in the matrix T. This proves the existence of T. 

We shall inductively prove the integrality of the matrix T using more explicit 
details of Eq. (13.111) . We begin with the equation for r = m — 1 and it gives 



7m t m— \"m— l,m— 2 "m,(m— 1) • 'm,m— 1 J- j v" - -'-^/ 

where we have used b m ^ m -\ = 1 for m > 1. Next consider, r = m — 2. This 
equation gives r mim _ 2 + r mim _i6 m _i im _ 3 = & m , m _ 2 which gives 

T~m,m—2 ^m^m—'l T~m,m— l"m— l,m— 3 > ^O.loj 

where we have used the fact that r miTO _i has been solved for and shown to be 
integral in the previous step. Note that this implies that T mm _ 2 is integral. 
Proceeding in this manner from r — (to — 1) to r — 1, we thus determine all the 
unknowns. A typical equation will take the form (reflecting the triangular nature 
of the equations) 



m—l 



b m ,m-r- £ r m ,A, 



' m,m—r 



(3.14) 



for r = 1, 2, . . . , (to — 1). We assume that T m , m _ r i is integral for all r' < r. Thus 
the right hand side is integral as it only contains integral terms. Hence r m ^ m - T is 
integral. This concludes the proof of integrality of the matrix T. 

□ 
We now state an unproven conjecture of Hanna and Meeussen. 

Conjecture 3.5 (Hanna-Meeussen) m\ r m+r+ i )m is a polynomial of degree m 
in r with integral polynomial coefficients. 

It is easy to show that r m+r+ i )m is a polynomial of degree (2m — 1) in r using the 
properties of the A-matrix. However, the above conjecture is stronger and seems 
to consistent with known data for m = 0, 1, . . . , 11. 

4 Practical Considerations 

This section provides details on the exact enumeration of higher-dimensional 
partitions as well as the triangles defined in this paper. With access to high- 
performance computing getting easier in recent times, this is indeed an additional 
computational aspect that can and must be added to the theoretical discussion 
of the previous section. We will first discuss the algorithms that we used and 
then discuss exact enumerations as we carried out. 
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4.1 Algorithms for higher-dimensional partitions 

There are two algorithms in the literature for computing higher- dimensional par- 
titions. The first one is due to Bratley and McKay (the BM algorithm) [7] and 
the second one is due to Knuth [8] - both are more than 40 years old reflecting 
the lack of progress in this area. Both are highly recursive and provide distinct 
ways of exactly enumerating higher dimensional partitions. 

The BM algorithm 

The partitions in any fixed dimension, say d, form a tree which we call the 
partition tree in (d+ l)-dimensiongj and denote by the symbol Td+i- Every node 
of the tree is the Ferrers diagram associated with a partition. The unique Ferrers 
diagram containing one point is the root node of the tree. New partitions can be 
formed by adding or deleting a point from the Ferrers diagram_|. Add a link to 
partitions connected this way. The depth of the tree is the number of points in 
the partition. 

The BM algorithm recursively traverses the tree up to some fixed depth, say 
n, such that each node is visited precisely once. The heart of the algorithm is 
the routine called part that takes three arguments and is recursively called in the 
algorithm. Every time a node is visited, the partition is stored in an array called 
current and presented to user. If one is interested in only counting the number 
of partitions of an integer in a given dimension, if the current partition has m 
points, increment a suitable counter, call it pdijn), by one. At the end of the 
program, the counter thus contains the number of partitions of all integers less 
than or equal to the depth of the traversed tree. 

The Knuth algorithm 

Let S m = N m denote the set of points in the totally positive orthant in a hyper 
cubic lattice. Let d m (k) denote the number of topological sequences with index k 
(see [6j|8] for definitions). Then a theorem due to Knuth |8j relates the numbers 
of topological sequences to numbers of partitions. To be precise, one has 

n 

p m (n) = ^2d m (k)p 1 (n-k) . (4.1) 

fc=0 

Since one-dimensional partitions are easily enumerated from the generating func- 
tion, it is simple to generate p m {n) given d m (k) for all k < n. Knuth provided 
an algorithm to generate and count all topological sequences - he illustrated this 



2 Recall that the Ferrers diagram for a d-dimensional partition is a set of points in d + 1 
dimensions. 

3 To avoid confusion, in this section alone, we shall refer to nodes of a partition as points in 
the Ferrers diagram. This is to avoid confusion with the node of the tree. 
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method by generating numbers for the numbers of solid partitions for integers 
< 28. Recently, a parallelized version of this algorithm was used by the author 
and other collaborators to enumerate solid partitions of integers < 68 [6]. 

Remark: An important aspect of the BM algorithm is that its memory usage 
is of the order of nd bytes, where d is the dimension and n is the maximum 
depth. This is vastly superior to the Knuth algorithm, where a similar problem 
needs memory of the order of n d ~ x bytes. However, when memory isn't an issue, 
our implementation of the Knuth algorithm typically takes less time than our 
implementation of the Bratley-McKay algorithm. 

The modified BM algorithm 

We begin with the observation that a suitably chosen sub-tree of the partition tree 
in r-dimensions, T r generates all partitions that contribute to the r-th column of 
the A- matrix i.e., o n>T .. The head node of this sub-tree is the Ferrers diagram \x T 
defined in Eq. ( I2.9J1 . The rest of the tree is generated by adding points to fi r . 
Let us denote this sub-tree by V r and the depth of this tree is clearly m where 
m — n — r — 1. 

The BM algorithm was designed to recursively traverse the partition tree vis- 
iting each node precisely once. The starting point of the algorithm is the root 
node whose Ferrers diagram consists of one point. Our idea is to change the ini- 
tial configuration in the BM algorithm to the Ferrers diagram, fi r and then call 
the recursive routine part with suitably chosen argumentsj. For this modification 
to work correctly, the program should traverse the sub-tree V r visiting each node 
precisely once to the chosen depth. This turned out to be easier as we exper- 
imentally observed that the sub-tree V r appeared naturally in the original BM 
algorithm for low values of r. We then checked that the modified BM algorithm 
correctly generated entries in the A-matrix for r < 10. However, we have not 
rigorously proved that this is indeed the case. 

Thus, once we have the modified BM algorithm correctly traversing the sub- 
tree V r , we can do the following: 

• Count the number of nodes at each depth - this gives the number a m+r+r ^ r . 

• At each node, numerically compute the reduced dimension, x of the Ferrers 
diagram. Then organizing the partitions by depth and r.d., we determine 
( r )cm,x- The binomial pre factor is present since all x < r will appear. This 
also implies that the algorithm is inefficient computationally for obtaining 
entries in the C-matrix. 



4 We have determined that the correct call is part [r + 2, 0, ( r ^ )). For comparison, the BM 
algorithm begins with the call part (1,0, 1). We thank Arun K. Jayaraman for implementing 
the BM algorithm as well as working out this modification. 
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A wish list of algorithms 

As we just mentioned, the current algorithm to enumerate entries in the C- 
matrix is computationally inefficient as we generate ( r ) partitions for each distinct 
contribution to c m>x . It is also inefficient because we need to compute x for every 
given partition. Can we create a more efficient algorithm? The problem is that 
we do not have an elegant characterization of sFD's with r.d. equal to x. This 
is in contrast to what happened with the A-matrix. In that case, we could show 
that any FD that has i.d. r necessarily contains the FD \x r . By using it as 
our initial configuration, we directly avoided configurations with smaller intrinsic 
dimension. For the C-matrix, we cannot avoid configurations that have smaller 
r.d. than the one of interest. 

We do not have any algorithms for the a and (3 matrices as well as the D/F 
matrices. So far these have been computed only indirectly after the A and C 
matrices have been computed. However, Proposition 12.161 might be a good start- 
ing point to coming up with an algorithm that directly enumerates entries in the 
F-matrix. 

4.2 Exact enumeration of higher-dimensional partitions 

In order to evaluate higher-dimensional partitions for integers < 25 and dimen- 
sions < 10, we chose to use the Knuth algorithm do carry out our computations. 
There was no serious memory issues for dimensions < 7 and the Knuth algorithm 
worked well. 

We needed to modify our computation when for dimensions 8, 9 and 10. The 
reduction in memory was done by counting topological sequences that fit into a 
box of size b. Then the memory requirement went down from n d ~ x to b d ~ x . For 
instance, when n = 20 and b = 10 (for d = 10), the memory usage went down by 
a factor of 2 9 and enabled us to keep our memory requirements in the 4 — 8 GB 
range as constrained by the IITM supercluster. However, some configurations are 
missed out as they do not fit into the box. Interestingly, one can show the error 
due to missed configurations is independent of box size when the index lies in the 
range [b + 1,26]. This makes it easy to estimate the errors by comparing with 
known results at smaller values of b and then slowly increasing the value of b. 
This method was used, for instance, to determine the ten-dimensional partitions 
of 20 - this was carried out by using a box of size 11 with errors determined up 
to k = b + 9. This was one of the more difficult computations as it took a several 
months of computer time to first estimate the errors and then carry out the final 
run in the box. Table [1] gives the results obtained used the Knuth algorithm for 
n < 23 and d < 10 and represent more than six months of computer time. 
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4.3 Exact enumeration of the A and C triangles 

The modified BM algorithm was used to generate the A and C matrix. The first 
eight rows of the C-matrix have been completely determined. Two additional 
rows were determined using additional information from the Z)-matrix. We obtain 

C m ,2m-2 = 6(2^-4)!! i 3 ^ ~ m ~ l ) 

c m , 2m - 3 = 6 gr 3 4)!! ( 2 ™ 4 - 6 ™ 3 + 3m 2 + 3m + 4) (4.2) 

C ™,2m-4 = 18 ( o 2 (^-6)!! ( 15m " ~ 75 ™ 4 + 95 ™ 3 + 21m ' + 88m + 42 ) 

C m ,2m-5 = 90 2 2 T-6)!! ( 258 ~ 167m ~ ^^ + lllm3 ~ 174m4 + llQm ^ ~ ^™ & + ^ 

This determines all entries in the A- matrix of the form a m+r+ i^ r for m — 0, ... 10 
for all values of r. We have determined the remaining entries for a n ^ r for n < 23 
by using the BM algorithm when necessary. The entry 023,11 was one of the 
longest runs and took about 880 hours of CPU time. Tables [2] and [3] provide our 
results. 

Using the (3 and a matrices as well as the Meeussen conjecture, we have also 
determined the 11-th row of the C-matrix. While none of these results were used 
in finally determining the entries in the A-matrix, there doesn't seem to be an 
inconsistency. This is only to be viewed as evidence for various conjectures. 

4.4 Extracting the elements of the other triangles 

All other triangles were obtained by using known numbers for the A and C 
matrices as we do not have an algorithm to enumerate them. The results for the 
D and the /9-matrix are presented in Tables |5] and [6] respectively. 

An improved implementation of the Bratley-McKay algorithm was provided 
to us recently by Prof. Bratley. This enabled us to enumerate a few more terms 
- in particular, we were able to enumerate rows 24 and 25 up to and including 
°25,i2- This enabled us to completely determine 25 rows of the F- matrix. This 
in turn determines all entries in 25 rows of the A-matrix and hence determines 
partitions of 25 in any dimension. It also provides a check on the 23 rows of the 
A-matrix which was independently determined. Table 0] provides our results. 



5 Concluding Remarks 

We have shown the existence of several structures that lead to simplifications 
in the exact enumerations of higher- dimensional partitions. The combinatorial 
interpretations that we have provided have enabled us to come up with an al- 
gorithms to evaluate the A and C matrices. A few lines of code in Mathemat- 
ica/Maple/Maxima/java can be used to store the A matrix and compute pd{n) 
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for n < 25 using the Binomial transform in real time. A working implementation 
of this is provided on the webpage: 

http : //www . physics . iitm . ac . in/~suresh/partitions . html . 

We will be adding these numbers to the OEIS as well as providing modules for 
SAGE/Mathematica/Maxima. 

It appears difficult to improve on our results which have determined all entries 
for the n = 25 row of the A-matrix. In fact, we have determined most of the 
entries for the n = 26 entry and hope to add this row in the future. Further 
additions to the A-matrix will require new and efficient algorithms to directly 
enumerate either the C or the F matrix. Another approach would be a naive 
parallelization of the BM algorithm. We hope to be able to eventually determine 
partitions of integers less than 30 in any dimension in the future. 
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A Ferrers Diagrams in a symmetric box 

Let us consider Ferrers diagrams of i.d. r that fit in a symmetric box of size b 
- points that lie within the box are such that all their coordinates take values 
in (0, 1, . . . , b — 1). Let us call them restricted Ferrers Diagrams. It is easy to 
see that under the action of S r that permutes the r-axes, FD's that fit in a box 
get mapped to FD's that also fit in the same box. Due to this property, we can 
construct analogs of the various triangles A/C/D/F for restricted FD's as well 
even though the total number of restricted FD's are finite. For instance, we have 



d+l 



rest 

Pd 



» = E( d r 1 )< ? ' (A ' 1} 
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where p T ^ st {n) is the number of FD's with n-nodes that fit in a symmetric box of 
size b d+1 and A rcst = (a^°f). Similarly, we can define C rcst . The analog of the A- 
matrix for restricted partitions not necessarily in a symmetric box has appeared 
in the work of Destainville et. al. [14j. 

Let us focus on partitions that fit into a symmetric box of size two and denote 
the corresponding triangles in obvious notation: 

A hox2 = (a b ° x2 ), C hox2 = (cjg*), D hox2 = « ox2 ) and F box2 = (/ b ° x2 ). 

We do not write out their relationships as they exactly mirror the corresponding 
formulae for unrestricted partitions. 

Definition A.l Let B be the set of strict Ferrers diagrams that fit in a symmetric 
box of size 2 . 

We say, somewhat loosely, that a strict skew FD, a of r.d. x is in B if the FD 
fi x U a G B. The only irreducible strict skew Ferrers diagrams at dimensions 2, 3 
and 4 in B are <T2, a?, and o^a/h respectively. It is also easy to see that V C B. 

The matrix C hox2 = (c^ 2 ) (for m, x > 0) has non-zero entries when m e 
[0, 2 x -x-l] with c£° x2 = 1. Further, c h Q °f = c h x °f = for x > and c^.^ = 1. 
Thus, it is a triangle. The maximum value of m, for fixed x, is obtained by 
considering the FD containing all nodes that are in the box. Such an FD has 2 X 
nodes in x-dimensions and thus the corresponding skew FD has 2 X — x — 1 nodes 
after deleting the nodes that lie in \i x . Below, we provide the first few rows of 
the matrix, C hox2 , for m e [1,6]. 



1 

3 3 



C b ° X2 = 1 16 30 15 (A.2) 

I n n n 1 1K iqc qqa qi k iak I V / 



1 15 135 330 315 105 
18 232 1581 4410 5880 3780 945 
. 13 355 4000 23709 71078 116550 107100 51975 10395 . 



Extending the ideas that were used in defining the D/F triangles, we look 
to count only those skew FD's that do not contain skew FD's in B are reducible 
components. Let the matrix C m>x denote this reduced C-matrix that counts strict 
skew FD's of r.d. x with m-nodes. Such FD's necessarily contain at least one 
node of type 2. Then one has the following relation that relates C to the C: 

x m / \ 
C m,x = / j / J II c m -p,x-y C P,y ' (A-o) 

y=0 p=0 ^' 

with co,o = 1 and co,a; — for x > 0. It is better to rewrite the above formula as 
follows: 



x—l m—1 



x—L m— i / \ 

c m,x + c m ,x "i / 4 / _, I I c m -p,x-y C P,V • \-™-^) 

,;=1 n=1 \y / 



C"m,x ^m,x i ^rn x 



y=\ p=l 
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The first term in the right hand side of the above equation is the contribution 
from skew FD's that do no contain any reducible components in B, the second 
term arise solely from terms that fit into a box of size 2. The last terms runs 
over terms that contain reducible components in B but do not fit into a box of 
size 2. The next proposition shows that C is a lower-triangular matrix with the 
m-th row containing m non-zero terms. 

Proposition A. 2 c m ^ x = when x > m or equivalently when the density p < 1. 

Proof: Since T> C B and all the irreducible strict skew FD's with density less 
than unity lie in D, the above Proposition follows from Proposition 12.131 □ 



C 



/oi \ 

(A.5) 



13 

17 16 

1 11 57 125 

1 18 135 602 1296 

1 26 293 1911 7980 16807 

\ 1 38 574 5242 31860 127977 262144 / 

\0 1 52 1089 12972 106505 619872 2411416 4782969/ 

\m— 1 





r-1 n ~ 1 /^\ 




Q"n,r " 


s =0 p=s+l \ ' 


box2 7 
u n-p+l,r-s JP,s 



We observe that c m , m = (m + 1) 

We can carry out a similar refinement for strict FD's that contribute to the 
A-triangle. One has 



(A.6) 



with fifi = 1 and / n .o = for n > 0. In order to interpret the first term, it is 
better to think of a n<T as the number of skew FD's obtained after removing the 
node at the origin of a strict FD. Then, the second term is the contribution from 
such skew FD's that do not contain reducible components that fit in a box of size 
two. A second equivalent definition in terms of m is as follows: 

r— 1 m / \ 
a m+r+l,r = Jm+r+l,r + a m _|_ r +i ir + /, , /, , I I a m— p+r— s+l,r— s Jp+8+1,8 ■ \™-') 

g=l p=0 ^ ' 

It is easy to see there is a bijective map that relates skew FD's that contribute 
to %a,x an d those that contribute to a m+x+ i jX . The bijection follows by observing 
that if a is a strict skew FD with m nodes and r.d. x, there is a unique FD (with 
i.d. and r.d. equal to x) obtained by adding the nodes in \i x . Thus, 



/, 



m+x+l,x 



It is easy to see using Proposition IA.2I that for / n>r = when r > n/2. We 
define the matrix F = (f n ,r) for n = 1,2,... and r = 0,1,2,.... Further, we 
observe that f2x+i,x — c x,x = {x+ l)*' 1 . Below we reproduce the first eleven rows 
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of the F-matrix. We reproduce the F-matrix alongside for comparison. The first 
instance where they differ is when n = 8 and r = 3 - this is precisely where the 
node (1, 1, 1) T that is not in V but present in B appears. As we go to higher 
values of n, an entry in F will be generically smaller than the corresponding entry 
in F. Further, entries in both matrices will agree when the density is in the range 
[1, 4/3) - this is because a node present in B but not in V first appears at density 
4/3. 



(« 


^ 




(* 


1 






1 


1 






1 


13 






13 


17 




, F = 


17 


1 11 16 






1 11 16 


1 18 57 






1 18 58 


1 26 135 125 






1 26 135 125 


V 1 38 293 602 
\0 1 52 574 1911 


1296/ 




V 1 38 293 618 

\ 1 52 574 1927 1296 





\ 



/ 



(A.8) 



The second row has only vanishing entries. That is because the only strict FD 
with two nodes fits in a box of size two. So the first non-vanishing contribution 
appears at n = 3, r = 1 if we ignore the n = 1, r = term that is more less part 
of the definition. 

We can now revisit the problem of enumerating partitions of n in any dimen- 
sion. We see that we need to enumerate the first n rows of the F matrix and A hox2 
in order to obtain row n of the A-matrix. However, from Eq. (I A. 6 1) we see that it 
is sufficient to determine only the first [n/2] elements of row n as that completely 
determines row n of F. However, this reduction is accompanied by the need to 
evaluate A hox2 which is yet another computation. Hence, we preferred to work 
with the F-matrix. However, one should be open to using the F-matrix if one 
has an algorithm to directly compute it. Then, the additional effort to compute 
A box2 might be worth it. 



A.l The box transform 

Define the following generating function for the A-matrix 

q m t r 



A{q,t) = ^^a m+r+ i ir 



m=0 r=0 



T\ 



(A.9) 



along with similar definitions for A hm2 (q,t) and F(q,t). Then, Eq. (1A.6I) implies 
that the generating functions have a simple relation. One has 



A(q,t) = A bo ^(q,t)xF(q,t). 



(A.10) 



It is due to this property that we refer to Eq. ( 1A.9J) as the box transform. Simi- 
larly, one defines 



C(q,t) = 5^5Z c "v 



q m t r 



m=0 r=0 



Tl 



(a.ii; 
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along with similar definitions for C hox2 (q, t) and C(q,t). Again, one has 

C(q,t) = C h °* 2 (q,t)xC(q,t). (A.12) 

There is an obvious extension to our considerations by replacing the symmetric 
box of size two by one of size b. Again, relations of the kind that we considered 
between FD's that fit in the box and those that don't appear. For instance, one 
has 

A(q,t) = A hoxb (q,t) x F(q,t) , (A.13) 

where A(q,t) is the generating function of FD's that don't fit into a box of size 
b and do not have reducible parts that fit into the box. 
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d = l 


d = 2 


d = 3 


d = A 


d = b 


d = 6 


d = l 


d = 8 


d = 9 


d= 10 
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10 


11 


3 


3 


6 


10 


15 


21 


28 


36 


45 


55 


66 


4 


5 


13 


26 


45 


71 


105 


148 


201 


265 


341 


5 


7 


24 


59 


120 


216 


357 


554 


819 


1165 


1606 


6 


11 


48 


140 


326 


657 


1197 


2024 


3231 


4927 


7238 


7 


15 


86 


307 


835 


1907 


3857 


7134 


12321 


20155 


31548 


8 


22 


160 


684 


2145 


5507 


12300 


24796 


46209 


80920 


134728 


9 


30 


282 


1464 


5345 


15522 


38430 


84625 


170370 


319555 


565983 


10 


42 


500 


3122 


13220 


43352 


118874 


285784 


621316 


1247780 


2350183 


11 


56 


859 


6500 


32068 


119140 


362670 


953430 


2240838 


4821737 


9661465 


12 


77 


1479 


13426 


76965 


323946 


1095430 


3151332 


8011584 


18478640 


39401792 


13 


101 


2485 


27248 


181975 


869476 


3271751 


10314257 


28395213 


70261505 


159527302 


14 


135 


4167 


54804 


425490 


2308071 


9673993 


33457972 


99845553 


265266530 


641733862 


15 


176 


6879 


108802 


982615 


6056581 


28310881 


107557792 


348333411 


994606250 


2565774277 


16 


231 


11297 


214071 


2245444 


15724170 


82033609 


342732670 


1205925033 


3704360354 


10198601886 


17 


297 


18334 


416849 


5077090 


40393693 


235359901 


1082509680 


4142850423 


13705110470 


40305279454 


18 


385 


29601 


805124 


11371250 


102736274 


668779076 


3389190112 


14122999548 


50367905030 


158376907546 


19 


490 


47330 


1541637 


25235790 


258790004 


1882412994 


10518508294 


47772540002 


183864216415 


618742851276 


20 


627 


75278 


2930329 


55536870 


645968054 


5249817573 


32361863632 


160336300356 


666612686420 


2403142436321 


21 


792 


118794 


5528733 


121250185 


1598460229 


14510628853 


98711666690 


533909133114 


2400146830007 




22 


1002 


186475 


10362312 


262769080 


3923114261 


39762851345 


298546248070 


1763901729589 


8581152930795 




23 


1255 


290783 


19295226 


565502405 


9554122089 


108058883583 


895425789360 









Table 1: d < 10- dimensional partitions of n < 23 as determined by direct enumeration using Knuth's algorithm, 
provides an independent cross-check of the entries in the first 11 columns of the A- matrix. 
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1 


40 
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54 
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3396 
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3514 
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1 


75 
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35115 


24318 
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55 


1 


13 


1 


99 


2185 


17910 


69580 


138155 


138075 


65997 


14667 


1495 


66 


14 


1 


133 


3765 


38942 


191795 


495870 


677663 


471276 


161202 


26875 


2156 


15 


1 


174 


6354 


82338 


505640 


1657975 


2978735 


2864408 


1424142 


360940 


46596 


16 


1 


229 


10607 


170265 


1285754 


5240090 


12016809 


15354492 


10604286 


3880561 


751696 


17 


1 


295 


17446 


345291 


3173220 


15821657 


45268685 


74497870 


68869266 


34954135 


9685709 


18 


1 


383 


28449 


689026 


7637795 


45999383 


161270025 


333494972 


400292769 


272579245 


104184949 


19 


1 


488 


45863 


1355253 


17996010 


129560563 


548523528 


1397398036 


2123894171 


1886698315 


965585764 


20 


1 


625 


73400 


2632975 


41631740 


355205608 


1794375520 


5541288850 


10446368715 


11819801575 


7897875909 


21 


1 


790 116421 


5058305 


94786545 


951526108 


5678296645 


20973892932 


48206965521 


68073453307 


58101011914 


22 


1 1000 183472 


9622420 212812255 2498219985 17463026868 


76290515426 210725428060 


364964576905 


390349624764 


23 


1 1253 287021 18139620 471921560 6444739208 52390397612 268136421612 879260678868 1840128105650 2425318710876 



Table 2: The first 12 columns and 23 rows of the triangle A. The other 11 columns can be obtained using the ten rows of 
the C-matrix given below. Thus, one can determine partitions of positive integers < 23 from it. 
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36285 
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5051970 
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945945 


135135 
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52 


1092 


13226 


114220 


751639 


3854487 


15231326 
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97613505 


150613155 


162889650 


120270150 


57702645 


16216200 


2027025 
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73 


1963 


30648 


332035 


2747799 


17918432 


92357844 


370929320 


1136808010 


2609559315 


4427605050 


5488733250 


4892112225 


3047969925 


1259458200 


310134825 


34459425 
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1 


97 3471 


67868 


910729 


9268382 


74767133 


483797592 


2498431224 


10155656364 


31998207087 


77214286182 


141528086700 


195617897475 


201837365730 


152796603960 


82323566325 


29876321475 


6547290750 



Table 3: The second triangle 
entries. 



the first ten rows and nineteen columns of the C-matrix. We have only shown non-zero 
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14200018 
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28205459 


65151254 


53672292 
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14617100 


86238209 


263040064 


359302890 
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34962755 
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975528302 


1899997612 


1309707840 
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4765327 


81792100 


711409264 


3398678150 


8749699709 


10128660960 


2357947691 


22 


1 


998 181102 


9157550 


187791450 


1948153500 


11278286646 


36739765288 


61114773760 


35600917115 


23 


1 1251 284021 17406714 


424233500 


5203415684 


35979941641 


144179174632 


318163092360 


314636749085 61917364224 


24 


1 1571 442713 32771292 


944990470 13605818265 111092074842 


536798419714 1499829016296 


2148096711540 1066426694784 


25 


1 1954 685443 61158328 2079070155 34930133300 333670251012 1915118952548 6574308285588 12551603978445 10672681371264 1792160394037 



Table 4: The F-matrix as determined using data up to 025,12 • 
only shown non-zero entries. 



This determines partitions of all integers < 25. We have 
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59024 
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1 38 


574 


5272 


33340 


154654 


509912 


1089488 


1158192 


378000 
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1 52 1089 12998 108465 


671389 


3123477 


10485214 


23226165 


28428750 14206500 1247400 


9 


1 73 1960 30336 321130 2551119 15580292 


72440912 


245511503 


561332710 762518790 501491760 102702600 


10 


1 97 3468 67430 891114 8811002 67908953 409620720 1895816757 6456110604 15166699372 22350118032 17852174340 5864859000 340540200 



Table 5: The first eleven rows of the D-matrix. We have only shown non-zero entries. 



CO 



Ay 







1 


2 


3 


4 


5 


6 




7 




8 


9 


10 





1 




































1 


1 




































2 


-1 




































3 


-5 


6 

































4 


1 


74 

































5 


41 


252 


-40 






























6 


31 


-540 


-676 






























7 


-461 


-6470 


1180 


656 



























8 


-895 


-4074 


69020 


10864 



























9 


6481 


138264 


403620 


-39016 


-9216 
























10 


22591 


376872 


-961240 


-1628984 


-191456 
























11 


-107029 


-2922930 


-21162456 


-3687040 


463680 


195840 





















12 


012.0 - 


-15867390 


-40350840 


168546560 


40336016 


7455104 





















13 


013 





013,1 


758778240 


1656046448 


110435472 


73922176 


-6297600 


















14 


014 





014,1 


014,2 


-1927766192 


-5730022032 


-552798336 


-382393600 


















15 


016 





015,1 


015,2 


015,3 


-44646818832 - 


-10585577760 


-7549384960 


278906880 













16 


016 





016,1 


016,2 


016,3 


016,4 


75450085920 


-14753227264 


22686050304 













17 


017 





017,1 


017,2 


017,3 


017,4 


017,5 


1603141023616 


607200778752 


-14729379840 








18 


018 





018,1 


018,2 


018,3 


018,4 


018,5 


018,6 


2727351931392 


-1449282760704 








19 


019 





019,1 


019,2 


019,3 


019,4 


019,5 


019,6 


01E 


,7 - 


-47662776674304 


873791815680 





20 


020 





020,1 


020,2 


020,3 


020,4 


020,5 


020,6 


02C 


,7 


020,8 


101710939668480 





21 


021 





021,1 


021,2 


021,3 


021,4 


021,5 


021,6 


021 


,7 


02 


,8 


021,9 - 


-58358690611200 



Table 6: The /3-triangle to the extent that we have determined it. The first column is consistent with the Meeussen 
conjecture. The (3 m+XtX for x G [0, m — 1] completely determine the degree 2m polynomial g m { r )- Thus, we have determined 
all polynomials for m < 11 albeit assuming the existence of the /3-matrix which is conjectural. The polynomials obtained 
this way agrees with the ones determined by the C/D/F matrices. 



m\x 


12 3 


4 5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 


17 


18 


19 20 



1 


1 

1 
































2 


3 


3 






























3 


1 


16 30 


15 




























4 


1 


15 135 


330 


315 


105 
























5 





18 232 


1581 


4410 


5880 


3780 


945 




















6 





13 355 


4000 


23709 


71078 


116550 


107100 


51975 


10395 
















7 





10 450 


8075 


78725 


431460 


1353240 


2552130 


2962575 


2079000 


810810 


135135 












8 





6 530 


14065 


204540 


1767045 


9207945 


29811330 


62179425 


85270185 


76621545 


43513470 


14189175 


2027025 








9 





4 580 22315 


456400 


5704580 


44793784 


225211165 


746795775 


1680747090 


2612970360 


2812925115 


2062160100 


983782800 


275675400 


34459425 




10 





1 611 


33177 


918981 


15738310 


174240318 


1268511894 


6207749790 


20975922462 


50107517460 


85928953110 


106306245045 


94166932860 


58305347100 


23983759800 


5892561675 654729075 



Table 7: The first eleven rows of the C box2 -triangle. We have only shown non-zero entries. 



CO 



n\r 


12 3 4 


5 6 


7 


8 


9 10 


11 


1 
2 


1 













3 















4 















5 















6 















7 















8 


1 












9 















10 


12 










11 


12 










12 


10 150 










13 





6 330 










14 





4 485 


2160 








15 





1 570 


7750 








16 





1 610 


17280 


36015 






17 





600 


30120 


185430 






18 





580 


45720 


574280 


688128 




19 





530 


63870 


1364195 


4727520 




20 





470 


85325 


2751875 


19192880 14880348 




21 





387 110625 


4994640 


59080000 130094748 




22 





310 140322 


8480885 


152220320 664737850 


360000000 


23 





215 174380 13808620 346973284 2557358244 


3873139200 


24 





155 212815 21879725 726316080 8167776498 24169328400 9646149645 



Table 8: The first 24 rows of the F box2 -triangle. We have only shown non-zero entries except for rows which have all zeros 
where we shown a zero in the first column. 



